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SECTION I 


HraRQDUCTION 


This report documents the development of a low-level approach system 
model, as one aspect of a terminal area system model, and includes an 
example of its application. The emphasis is on point of view, technique, 
and simplicity. As such, the report takes on the characteristics of a 
tutorial presentation. 

A. RJRPOSE 

The primary purpose of a terminal area system model is to establish a 
structure containing the system elements, coramand inputs, disturbances, and 
their interactions in an analytical framework so that the relative effects 
of changes in the various system elements on precision of control, pilot/ 
copilot workloads, and available margins of safety throughout the terminal 
operations envelope can be estimated. It is intended that any such model 
will be used to provide insight for the design and integration of suitable 
autopilot, display, and navigation elements; and to assess the interaction 
of such elements with the pilot /copilot, with emphasis on the automatic/ 
himan interfaces. 

The model presented herein is an analytical one, and thus has numerical 
measures of performance as outputs. These outputs include the performance 
of both human and inanimate components, as well as of the system as a whole 
Thus, the model is expected to be extremely useful in: 

1 . Identifying those system areas which offer the 
largest immediate possibilities for improvements 
in safety of operations (and those that offer very 
little possibility for improvement). 

2. Quantifying the prediction of relative performance 
and safety (or success) margins for competing system 
(and subsystem) alternatives. 

3 . Identifying needed research to (a) fulfill a system 
need, or (b) improve the accuracy of the model. 
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h . Providing a long-range potential^ >^hen the model 
is verified;, for computing absolute pei-formance 
and safety limits (as opposed to relative levels 
among competing systems) to serve as a guide to 
the specification of subsystem requirements and 
as a means for estimating operational statistics. 

B. POINT OF VIEW 

A logical starting point for presenting the model is to define the 
desired outputs. Because the approach and landing phases of flight are 
the most critical, these phases were selected for analysis. This decision 
led to the desired outputs being the ’^Basic Outcomes” shown in Table I. 

In the table the Basic Outcomes are listed together with their "Associated 
Performance Measui^es" and "Performance Metrics" (which are computed from 
the dynamic portions of the system model). A key point in the model is 
that "Outcome Probabilities" can be computed by combining the critical 
limits on the performance measures with the values of the performance 
metrics. Although these probabilities will be only rough approximations 
for the first-cut models they will be useful as gross indices and will 
allow valid comparisons among alternative systems. 

As noted in Table accidents are very unli>:ely events^ and therefore 
all accident probabilities are extremely small. Consequently;, from a 
practical standpoint;, it is almost a necessity to consider the probability 
of a missed approach as a primary measure of system adequacy. This point 
of view is supported by existing information indicating that^ at the 
Category II-B level (see Table II);, the missed approach rate can be as 
high as 40 percent (Ref. 1 ). In this report we have adopted the FAA 
position concerning missed approaches (Ref, 11). This is that an 
approach can be continued below the decision height* (assimiing the 
pilot can see to land) only if the airplane is within 12 ft vertically 
from the center of the glide slope beam^ within 72 ft laterally from 
the center of the localizer beam;, and within 5 kts of the nominal approach 


*See Table II for a definition of decision height. 
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CONSIDERATIONS FOR APPROACH AND LANDING OUTCOMES 
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TABLE II 


ICAO LOW-VIS IB ILITY- LANDING HE CATEGORIES 


CATEGORY 

RUNWAY VISUAL 
RANGE (RVR) 
ft 

DECISION 
HEIGHT (DH)* 
ft 

I 

2,400 

200 

II-A 

1,600 

150 

II-B 

1,200 

100^ 

III-A 

TOO 

— 

III-B 

150 

— 

III-C 

zero 

— 


A height (above runway elevation) belov/ v/hich a pilot must 
not descend if he has not obtained adequate visual references 
to land; i.e*^ he must execute a missed approach at this point 
if he does not have adequate references to land by visual means. 

^Sometimes called CAT II. 

speed. These limits constitute an effective "window" at the decision 
height. If the airplane is not within the window^ then a missed approacl^ 
is mandatory. For practical operations^ making the Category II-B window 
a very large proportion of the time^ say 9v“t>9 percent, may be the most 
difficult part of the approach and landing sequence. However, once this 
window is made, only the flare and removal of crab or wing-down are 
required to get on the ground with reasonable touchdown conditions . 
Therefore, the occurrence of short landings arid off-runway landings 
should decrease when using the Category II window (while most of the 
other outcome probabilities would be expectea to be fairly similar to 
current operational experience). In any event, for a first-cut evalua- 
tion, we will direct most of our calculations to considerations of suc- 
cessful approaches and missed approaches. This will in no way detract 
from the generality of the technique, but will greatly simplify the 
mathematics involved. Hence, a much more straight-forward presentation 
is possible. 
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C . OVERVIEW OF AMLYSIS 


To estimate the performance metrics and outcome probabilities 
requires a dynamic model and probability performance trees (among 
other things). An overview of the various analysis steps is pro- 
vided in Fig. 1 . The dynamic model is made up of blocks 1 to 4 . 

The results of these first blocks are^ in essence^ feedback control 
systems which satisfy the guidance and control requirements for 
approach and landing. The airplane is the controlled element and 
the active elements of the feedback controller may be either the 
pilots autopilot^ or a split-axis pilot /autopilot combination. In 
addition, the systems contain passive or monitoring elements: the 

copilot, for the manual control situation, or the pilot and copilot 
for the automatic condition. The system mechanizations (Block 2) 
can be made appropriate for: 

• Fully automatic approach and landing 

• Flight director plus pilot/vehicle approach 
and landing 

• VFR pilot /vehicle approach 

The first two system mechanization possibilities are representative of 
advanced low- visibility approach and landing systems. The last one is 
still semiconjectural in that a scanning multiloop pilot model has not 
yet been validated for such a complex situation. 

Exercise of the dynamic model with the inputs and disturbances 
selected in Block 4 provides the performance metrics of Table I. Con- 
ventionally, this block, noted as 7 on Fig. 1, would use ordinary air- 
craft motions and kinematic quantities as the dependent variables. 
However, in some cases we have found it possible to select combina- 
tions of these variables as composite "state" variables which are 
more directly related to the basic outcomes (such as runway overruns) 
than are the standard aircrajfb motion quantities . Consequently these 
state variables are used, where appropriate, as part of the basic dynamic 
model to ease the transition between performance metrics and the ultimate 
outcome probabilities. With the composite state variables included, the 
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Figiire 1 , Overview of Analysis Steps Leading to Evaluation 
of Category II Approach and I.anding Systems 












"extended" dynamic model then comprises Blocks 1 through 6, with Block 7 
bej.ng its actual exercise for a particular system and set of conditions. 

With the performance metrics available, and critical limits and 
approach outcome probabilities defined (Blocks 8 and 9 in Fig. 1), the 
numerical values for the approach outcome probabilities can be computed 
(Block 1 O) . The result is then a set of performance metrics and outcome 
probabilities for a specific system mechanization subjected to a given 
input and disturbance environment. Direct comparisons of automatic sys~ 
terns can be made simply by going through these analysis steps for each 
competing system mechanization with a standardized set of inputs . It is 
noted that the role of incorrect pilot decisions for the execution of 
missed approaches is not considered here. Such a refinement is left to 
future versions of the model. 

In this report the intention is to develop and describe the model as 
outlined above . 

D. OUTLINE OF THE REPORT 

Using the framework from the above overview, Sections II through V 
present the information called for in the various blocks of Fig. 1 . 
Section VI then follows with a summary and some conclusions. Peripheral 
material is found in the appendices and in Ref. 12. 
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SECTION II 


DESCRIPTION OF FUNCTIONS TO BE PERFORMED DURING APPROACH 

During the approach phase of flight there are a large number of tasks 
to be performed (and monitored) by the flight crew. These tasks include 
going through checklists^ tuning radios , carrying on communications with 
controllers on the ground^ navigating per ATC clearances^ etc., as well 
as flying the airplane. A list of these tasks, when they should be per- 
formed, and the associated numerical parameters are given in Fig. 2. In 
addition to this overall view of the approach situation, Fig. 3 has a 
detailed view of the situation in the vicinity of the airport. 

For our purposes, the tasks of primary interest are glide slope and 
localizer tracking from the outer marker to the decision height. These 
tasks can be performed automatically (with an autopilot), manually (with 
a flight director), or via a combination of automatic and manual control. 

In this report we will consider the functions which must be accomplished 
by any form of control, and shall illustrate their actual performance both 
with fully automatic equipment and a manually controlled flight director. 
Modifications made to the analytical details to make them suitable for 
manual control are described in Appendix B.* In t-he initial system analy- 
ses no autothrottle or speed control is considered. Thus the basic system 
is evaluated first, and the effects of system perturbations can be assessed 
later. A detailed discussion of the functional requ.iremerits of the control 
system and the means for satisfying the requiremt ‘ut ;; is presented in Ref. 12- 

Although only the primary tasks of glide slops and localizer tracking 
are considered herein, such secondary tasks as system monitoring and fault 
detection are recognized as being equally important for achieving a success- 
ful approach. However, analysis of these aspects of system operation is 
beyond the scope of this initial system modelling,. 


* Briefly, the major modifications required to cover manual control are 
the addition of pilot lags in responding to displayed errors, and the intro- 
duction of pilot -generated remnant into the elevar.or, aileron, etc., com- 
mands. The amount of remnant is determined by t:ie vii splays used and the 
associated scanning required to close the various loops . 
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Figure 2, TasI^s Performed during Approach and Landing 




150 kts 



Figure 3 . Typical numbers Associated with an ILS Approach and handing 


SECTION III 


SYSTEM INPUTS 

For the purposes of evaluating and comparing systems it is necessary 
to subject them to an appropriate set of inputs. To cover a variety of 
situations^ three types of input disturbances which are commonly encountered 
during an approach were considered. 

• Random gusts 

• Deterministic gusts 

• ILS beam noise 

They shall be described in this same order below^ and then summarized in 
a table at the end of this section. 

A- RANDOM GUSTS 

The random gusts considered here are defined by components along all 
three axes. Since analysis of measured gust data has revealed nearly 
Gaussian distributions^ the assumption of zero-mean Gaussian distributions 
for random gusts appears to be justified. 

Considerable gust data from numerous sources have been integrated to 
produce a gust model appropriate for design analysis purposes. The 
pertinent aspects of this model are given in Ref . 2 . The Dryden form 
of the model was selected for ease of computation. The basic form of the 
model is defined by three power spectral densities which relate normalized 
gust intensities to spatial "frequency." In terms of temporal frequency 
these spectra are given by 



n 



where Ir^, In,; 


scale lengths fft) 


Uq " aircraft mean speed wit/; respect to the 
air mass ( ft / sec ) 

fjo = frequency (rad/sec) 

a - standard deviation 

An additional spectrum of interest is derivable from these basic spectra^ 
and is given by 


$Pg{-) 



(U 


where b = wing span (ft). 

Conforming to the stip\ilations given in Rel . the random turbulent 
velocity components have been assumed to be uncorrelated. While this 
assumption is fiilly justified for clear air turbulence at high altitudes^ 
it is not strictly true at low altitudes because ol' prevailing anisotropy 
in the boundary layer. Ho\\»ever, it is also pointed out in Ref. 2 that at 
low altitudes the cross -correlations between the gust components are weak, 
and may therefore be disregarded. 

In viev/ of the above ^ it is assumed that the tbiree velocity components 
Ug, Vg, and v/g, are mutually uncorrelated^ so that analysis can be carried 
out using the three components of the gust model separately. 

The procedure for evaluating the spectra is outlined in the following. 
It is based on extensive data fitting and adjusting of the Dryden "scales" 
to m^ake all three scales (L^, Ly, and Ir,;) equal at an altitude of 1,750 ft 
The resulting scale lengths for clear air turbulence for the 
above spectral forms are 

Below 1,750 ft: - h(ft) ( 5 ) 

Lu = Lv = (6) 
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The variation of Lq and Ly a,t low eultitudes according to the one-third 
power of altitude above ground level is simply a mechanism that forces the 
scales of the two horizontal gust components to be larger than the vertical 
scale. (Although these formulae produce correct trends, there is little 
data available that can be used to substantiate the h ^ relationship in 
SJld. Lry • ) 

This gust model represents stationary, randan gusts as stationary 
processes- However, the gusts actually encountered by a descending air- 
plane are nonstationary due to the altitude dependence of the gust charac- 
teristics. Further, because the breai: points in the analytic expressions 
for the gust power spectra are nonlinear functions of altitude, the spectral 
characteristics of the gusts encountered by a descending airplane are non- 
linear functions of time. However, in spite of these complicating factors, 
reasonably accurate performance calculations can still be made quite simply* 
This is explained in the following paragraphs. 

It is a straightforward calculation to ccampute the rms output of a 
linear system to a stationary random input- Figure U illustrates this- 
Implicit in these calculations is the assumption that the system has been 
operating for a sufficiently long (theoretically infinite) time so that all 
transients have died out. However, for practical purposes all that is really 
required is that "steady- state’* conditions have been reached, i-e., that 
the system has been operating on the stationary input long enough for the 
output to become approximately stationary- 


2 

-00 

•'O 


Figure h. Computation of RMS Output from Input Spectrum 
and System Transfer Function 


Input Output 

Spectrum Spectrum 



Linear 

System 
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The time required for this to occur depends on the system dynamics. For 

our model airplajie plus controller this time is relatively short of the 

order of 5 to 10 sec (based on the settling time for step inputs). As an 
analogous example, Fig. 5 shows a plot from Ref. R of versus time for 
an f 4D-1 tracking a visual glide slope beam during a constant speed 
approach to an aircraft carrier. For this plot, random beam motion was 
initiated at t = 0 sec. The plot shows that 8o percent of the steady- 
state airplane is reached in 5 sec, and 90 percent in 7-1/2 sec. 



Figure p . Aircraft Altitude Excursions 
Resulting from Pilot Tracking the Beam 

Thus^ althou^ the mathematics requires a stationary input to be 
applied to a system for an infinite time, only the most recent 10 sec 
(or so) of input has any appreciable effect on the airplane’s current 
condition. This observation has profound consequences^ the most 
pertinent of which is that the complete non stationary gust input to the 
descending airplane is not required to estimate -fne airplane dispersions 
at the Category II decision height of 100 ft latitude. Instead, only 
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the gust characteristics corresponding to the last 5- TO sec preceding 
the altitude of interest need he considered. Because the airplane sink 
rate is about 10 ft/sec^ ^ to ^0 sec corresponds to about 50 to 100 ft 
altitude change. For such a small altitude change the altitude-dependent 
gust parameters do not change very much during the time interval of 
interest. Therefore^ the required "short-time” gust model is essentially 
a stationary input which is closely approximated by the gust model 
components evaluated at an altitude of 100 ft. This yields 


Iv = 100 (ft) (7) 

1^ = = U5(iOO)h5 = (ft) (8) 


The gust intensities along the various axes for Dryden form spectra 
are related by the expression 

^2 2 2 




Py 

(9) 

ThuS;, the ratios 

of 

and a., to a,, will be 


% 

. hi 

^Wg 

= = 2 59 

1 If 100 

(10) 


The probability of occurrence of the gust intensities a, at various 

Wg 

altitudes is represented by P(a^ defined as the exceedance probability 

o 

p(%) = PiP(%) (11) 


where 


P-j is the probability of occurrence of 
turbulence and 

A 

P( ) is the probability of equalling 
a^given magnitude of Oy once clear 
is encoxmtered . ^ 


clear air 


or exceeding 
air turbulence 
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The probability of occurrence of clear air turbulence^ , at various 
altitudes is defined by the curve in Fig. 6 derived from various measurements. 
According to Fig. 6, P-| - 0 .8 at 1 00 ft. 


The probability P( 
depicted by the curve in Fig 


of equalling or exceeding a given a is 

W rr 


B. DETERMHTISTIC GUBTS 


1 . Wind Shears 

A 1 though horizontal wind shear near the ground is a relatively 
cammon phenomenon^ it is still not fully understood. This is due^ 
in part^ to the paucity of shear measurements. As a result^ such 
things as probability of encounter^ and distribution of shear magnitudes 
at various altitudes are not presently known. Another reason for the 
lack of a better understanding of wind shears is that shear can result 
from a number of different causes. Probably the best londerstood cause 
of horizontal wind shear is the boundary layer effect of the ground on 
a moving air mass. Thus^ the air closest to the ground moves slower than 
the air higher up. This primary effect then leads to a secondary effect 
because of the so-called Coriolis forces. The net result is that the wind 
shifts in direction in addition to decreasing in magnitude. This is a 
familiar situation to student pilots who are taught that the wind direc- 
tion will shift about deg (counterclockwise) and the magnitude will 
drop by almost fO percent during the final two tliousand feet descent to 
a landing (Ref. f ). The explanation for this is that the wind tends to 
align itself with the pressure gradient (from high to low} near the ground^ 
and to align itself with the Coriolis-produced "cyclonic” swirls (that are 
perpendicular to the pressure gradient) at higher altitudes. Figure 8 
depicts this situation simply. 

A sample of British data (Ref. t) for shear measurements is presented 
next to indicate the measured probability distribution of shears over rela- 
tively flat terrain. Figure 9 shows a histogram of the average shear over 
a height range of If meters down to 1 meter. Superimposed on the histogram 
is a Gaussian probability curve having a mean and rras deviation from the 






Probability of Equaling or Ex 
Once Turbulence has 



Figure 7. Hxceedance Probability for 





Figure 8 . Idealization of Wind Directions Near the Ground 
and at High Altitude (in the Northern Hemisphere) 

mean equal to that computed from the data comprising the histogram. A 
similar figure for average shear over an altitude range of 75 meters down 
to 50 meters is given in Fig. 10. Figures 11 and 12 then follow with 
cumulative probability plots to test the data for being Gaussian. As 
seen in these plots ^ the data appears to be reasonably Gaussian over the 
ranges measured. 

It is pertinent here to make a few comments regarding wind shear 
measurements. Tower data (giving simultaneous wind speed measurements 
at several altitudes) has consistently led to smaller shear values than 
is computed from instrumentation aboard a descending airplane. Although 
a purist may regard the tower data as the more appropriate method for 
obtaining accurate shears^ it is really the effect on the airplane that 
is of interest. That is^ the quantity actually desired is the rate of 
change of wind along the airplane ^s flight path (and not the wind gradient 
measured vertically over a single point on the ground at a given instant 
in tine). Thus^ the technical problems of measuring "accurate" wind shears 
is confounded by a semantic problem as well. Strictly speaking^ a wind 
shear is the instantaneous vertical gradient of horizontal wind. But, in 
addition to this type of wind variation, an airplane may experience a 
change in horizontal wind due to a longitudinal gradient, or even due to 
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Figure 10. Probability Distribution of Measured Wind Shears fron Ref. 

^“50 '■ h -■ Y'; meters^ 




Cumulative Probability of Wind Shear (%) 



Figure 12. Cumulative Probability of Measured V/ind Shears 
from Ref. ^ (30 £ h < 7^ meters) 
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a wind fluctuation as a function of time (’which may have no spatial gra- 
dient)^ or to any combination of these possibilities. When thought of in 
this light it seems reasonable that the effect of a true vertical gradient 
may actually be minor compared to the other effects. However, whether it 
is or not is really not of concern here. V/hat is important is to recognize 
that tower data is not truly representative of the magnitude of horizontal 
wind changes that are commonly experienced by an airplane during final 
approach. 

It may be useful at this time to present some direct quotations from 
Ref. 6 pertaining to Lear Siegler, Inc., experience with automatic landings 
in the Caravelle. In particular, the following concerns the occurrence of 
large wind shears near the ground. 


'^The experience obtained during the LSi/SUD program 
indicates that this gradient in practice is much more 
severe than generally accepted, and it also appears that 
this gradient becomes more severe as the altitude decreases, 
due to the effects of the ground on the air mass. Fore and 

aft wind shears of 50 knots per 100 feet of altitude, lasting 

for periods of eight seconds, have actually been recorded by 
rather complete instrumentation on at least three of the 
automatic landings made in Toulouse. The terrain at Toulouse 
is relatively level, and would not be considered conducive to 
causing such wind shears. The weather conditions at the time 
these occurrences were recorded did not appear to be abnormal." 

"Experience on the LSi/SUD program has also indicated that 
wind gradients approaching the same magnitude appear in the 
lateral case. Crab angles of 15 degrees have been experienced 
at 150 to 200 feet of altitude with the touchdown occurring with 

a zero crab angle. At the approach speeds of the Caravelle, this 

is equivalent to a cross-wind gradient of 16 knots per 100 feet, 
which is far in excess of four knots per 100 feet," 

"The wind gradients actually experienced in the LSi/SUD program 
disagree with the information presented in RTCA SC -79^ which is 
the source of the presently accepted four knots per 100 feet. 

However, the larger figures obtained on the LSi/SUD program are 
not the result of an isolated occurrence; also, they are well 
documented, and as such should be considered valid," 

Taking into consideration all of the above information on shear, 
"standard" shear inputs were selected for the longitudinal and lateral 
situations. For our calculations, wind shear is simulated by introducing 
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an altitude dependent steady wind — for both the headwind and crosswind 
components. For the headwind component the variation with altitude 
starts at 200 ft and is linear down to an altitude of 100 ft. Changes 
in the magnitude of the headwind shear occur at 100 ft and 50 ft, as 
shown in Fig. 15. The crosswind component also starts at 200 ft and has a 
linear variation with altitude down to the ground, as also shown in 
Fig, Although the magnitudes shown may be considered somewhat 

arbitrary, they are representative of actual measured shears, and are 
consistent with current thinking in the industry concerning autopilot 
requirements for Category III conditions. 

Having presented a model for random gusts, then a short discussion 
on shear, followed by some measured shear data, and, finally, the comments 
of an autopilot designer and test pilot, it is now of interest to present 
a casual consensus of the aircraft industry's unofficial thinking on over- 
all environmental condition limits for automatic approach and landing sys- 
tems. Table III contains such a consensus. It is noted that "patches" of 
turbulence and discrete gusts (such as steps and 1 — cosine pulses) are not 
commonly included in the environmental conditions pertinent to approach and 
landing. 


TABLE III 

CONSENSUS OF ENVIRONMEINTAL CONDITION LIMITS 
FOR AUTOMATIC APPROACH AND LANDING 


PARAMETER 

LIMITS 

Steady Wind 

25 kt headwind 
10 kt tailwind 
15 kt crosswind 

Turbulence 

Moderate: 6-1/2 fb/sec — 9 ft/sec rms Wg 

1 

Shear ' 

1 

Ug: -4 kts/100 ft TOO ft < h < 200 ft 

-8 kts/100 ft 50 ft < h < 100 ft 

-25 kts/lOO ft 0 < h < 50 ft 

V : 15 kts/100 ft 0 < h < 100 ft 
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2. Step Gusts and 1-Cos Gusts 


In addition to random gusts and shears, steps and 1 — cosine gusts 
should also be considered. The 1 - cosine gusts are included to represent 
the large discrete gust pulses that occasionally occur in turbulent air, 
but that are not appropriately included in the random gust or wind shear 
model inputs. The same magnitude was selected for all three gust components 
(ug, Vg, Wg), namely, 15 kt at the peak. This corresponds to the U. S. 
Weather Bureau* s definition of moderate turbiolence (where peak gust magni- 
tude is between 20 ft/sec and 35 ft/sec). The duration of the 1 — cosine 
pulse is 2.5 sec, which corresponds roughly to the time to travel 
25 chord lengths (Ref. 7 ). 

Step gusts are included to represent the occasional very long lasting 
gusts . Because the longer lasting vertical gusts near the ground are not 
as large as the horizontal ones, and because the longer lasting horizontal 
gusts are caused primarily by changes in the "total" wind, the magnitudes 
selected for the step gusts are different for all three components. The 
Ug and Vg steps were chosen to be 12 kts and 7 kts, respectively — which 
corresponds closely to a ^0 percent decrease in the total wind magnitude 
for the case of a steady wind having headwind and crosswind components 
as given in Table III. The 7 kt crosswind component also agrees with the 
crosswind gust environment specified for the C-lkl in Ref. 8 , although 
for the C-l4l the requirement was for a steady 25 kt crosswind — gusting 
to 32 kts. The 5 kt w^ step was selected primarily on the basis of Ref. 9 
which suggested it as a "realistic condition." Note that this makes the 
w component the smallest of the three, per the earlier comment regarding 

o 

the relation between long lasting horizontal and vertical gusts. All of 
the selected inputs are summarized in Table iv at the end of this section. 

C . IIS BEAM NOISE 

The effects of beam noise, for both the localizer and glide slope 
signals, can be determined by exciting the system with beam noise 
inputs . A localizer noise input was obtained from average power 
spectral density plots of beam noise for localizers at several air- 
ports . Plots for conventional and directional localizer noise are 
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<t>(u;)(/iaVrad/sec) <J>(a»)(/i.aVrad/sec) 


repeated here as Figs. 14 and 15 for easy reference. These spectral den- 
sity plots were scaled so that their integrals give the mean- s quared value^ 
i.e,^ 



Using Eq. 12 and the plots in Figs. l4 and V; leads to the result 
that the noise for the ''average’' localizer is of the order of 5 ::a rms . 
This exceeds the Category II requirement of 2f; [.la rms^ and therefore 
cannot be used directly . However^ for model purposes a localizer noise 
power spectrum was obtained by fitting a simple analpd:ical expression 
to the shape of the directional localizer^ and then scaling down the 
rms level to 2.5 pa (as shown by the dashed line in Fig. 15 ). 

We were unable to find an existing power spectral density plot for 
glide slope beam noise. Therefore^ we obtained some glide slope beam 



Figure l4. Avora, ;e Conventional localizer Power Cpectral Density 



Figure 15 . Average Directional Focalizer Povrer Cpectral Density 
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data from the FAA and analyzed it ourselves to obtain a PSD plot. The 
glide slope noise data came fran measurements at LaGuardia made by the FAA 
in January of 1967. This particular beam just meets Category II criteria 
frcxn the outer marker down to the middle marker, and therefore the PSD plot 
pertinent to this segment of the beam is the one used here. The PSD plot 
is given in Fig. l6 along with a data fit by a simple analytical expression. 



Figure 16, Power Spectral Density Plot of Glide Slcpe Beam Noise 
at LaGuardia ( from Outer Marker to Middle Marker at Low Tide ) — 

RMS Level is 10 ua 

It is noted that a short-time- stationarity argument similar to that 
used earlier for random gusts can be applied to beam noise. That is, 
even though the noise properties may be rajige-dependent, they are 
relatively stationary over intervals of the order of airplane/system 
settling times, and can therefore be treated as stationary inputs. 

It is pertinent here to point out that criteria for maximum allowable 
beam noise are currently specified only in terms of mean beam angle, maxi- 
mum bend amplitude, and rms deviation frati the mean. However, for simula- 
tion piirposes it is necessary to know also the amplitude and frequency 
distributions of the noise. (For example, a beam may have a Gaussian 
amplitude probability distribution and a band- limited white-noise frequency 
distribution. ) 
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TABLE IV 
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BEAM NOISE 





SECTION IV 


DERIVATION OF APPROACH OUTCOME PROBABILITIES 


One of the key aspects in the presentation of a terminal area model 
is a definition of the model outputs. This was accomplished in Table I 
(repeated here for ease of reference) which lists the "basic outcomes" for 
an approach^ and their associated measures^ metrics^ and critical limits. 
The metrics are computed from the dynamic portion of the aircraft/approach- 
system model. Then^ by comparing the metrics with the appropriate critical 
limits^ figures of merit in the form of outcome probabilities are computed. 
The various outcome probabilities are a basis for comparing competing sys- 
tems^ as well as a means for evaluating the effects of inputs and/or 
changes in system components . 

Figure 1 7 shows an approach outcome tree which indicates how aircraft 
dispersions canlead to the various approach outcomes listed in Table I. 
Although this figure is somewhat oversimplified^ it does indicate concep- 
tually how a model can be constructed to give approach outcomes as outputs. 
Before proceeding to a more detailed breakdown of the models a digression 
to clarify the definitions of the accidents listed in Table I is pertinent. 

A. DETAILED DISCUSSION OF TYPES OF APPROACH OUTCCMES 

During an approach more than one of the basic outcomes listed in 
Table I may occur. For example, a hard landing and running off the runway 
could both occur on the same landing. However, for such a situation one 
might still consider this to be a single accident because it is not sur- 
prising to find that an airplane ran off the runway after a particularly 
hard landing (where the landing gear failed). Thus the "critical" outcome 
can be considered to be a hard landing (for this case) even though another 
of the basic outcomes also actually occurred. Another example of a com- 
bination of basic outcomes would be landing short and having an excessive 
crab angle (enough to break the landing gear). Clearly, the fact that an 
excessive crab angle was present would go unnoticed if an airplane touched 
down on somebody's garage roof a mile short of the threshold. Thus the 
critical outcome for a combination of these two basic outcomes is the short 
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CONSIDERATIONS FOR APPROACH AND LANDING OUTCCMES 
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Figure 17, Simplified Version of Approach Outcome "Tree"’ 
(General Breakdovm of Possible Approach Consequences) 





TABLE V 


DEFINITIONS OF "CRITICAL" APPROACH OUTCavlES* 


BASIC OUTCOMES 
(AND ABBREVIATIONS) 

CRITICAL OUTCa.ffiS RESULTING 
FROM ANY TWO BASIC OUTCOMES ^ 

Successful Landing 

— 

OK 

OK 









Successful Missed Approach 


— 

MA 








Short Landing 

SL 

— 

— 

SL 







Hard Landing 

HL 

■ 

■ 

SL 







Overrun Runway During Rollout 

OR 

I 

■ 

D 







Land Off Side of Runway 

LO 

— 

■ 

■ 







Drag a Wing Tip or Engine 
Pod During Landing 

DW 

— 

— 

SL 

HL 

DW 

+ 

OR 


DW 1 





■ 

■ 

SL 

i-IL 

V 


DW 

+ 

V 

V 


Run Off Side of Runway 
Dur'^ng Rollout 

RO 

— 

— 

SL 

rH 

RO 


DW 

+ 

RO 

V 

RO 



OK 

MA 

SL 

HL 

OR 

LO 

DW 

V 

RO 






:asic 

OUTCOMES 





*For situations where an approach results in more than one of the 
basic outcomes. 

tA basic outcome from the left column paired \^lth a basic outcome at 
the bottom gives the critical outcome shown in the box for that row 
and column. 














landing. On the other hand^ sane combinations of basic outcomes deserve 
joint consideration if they occur on the same approach. Table V presents 
the critical outcomes for all possible pairs of basic outcomes. A list 
of the different critical approach outcomes from Table V is then given in 
Table VI.* 


At this point it is convenient to 
neglect the combination of dragging a 
wing tip and overrunning the runway. 
This is not unreasonable because over- 
running a runway is usually so much 
more serious than dragging a wing tip 
that if the two occurred on the same 
landing the dragged wing tip could 
easily be lost in the confusion. The 
convenience arises because the out- 
comes can now be separated into longi- 
tudinal and lateral situations^ as 
shovm in Table VII. 


TABLE VI 

LIST OF POSSIBLE CRITICAL 
APPROACH OUTCOMES 


OK 

MA 

SL 

HL 

OR 

OR +DW 
LO 
DW 

V 

RO 

DV/ + V 
DW + R0 


Before getting back to the detailed breakdown of the models it is 
pertinent to present some of the results of our investigation of accident 
and incident statistics for several recent years, (A detailed presenta- 
tion of accident and incident statistics is given in Appendix A.) One of 
the purposes of the investigation was to determine the relative likelihood 
of the various outcomes considered in Table I, In this way we could learn 
which outcomes were most important and which ones were of least importance. 
Table VIII shows those outcomes from Table I for which data is available^ 
and their relative likelihoods (based on the landing accidents and incident 
for 196U through 1966 for U. S. air carriers). 


*Our selection of critical outcomes is done to make our outcome 
definitions consistent with those used in publishing accident and incident 
data. 
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TABLE VII 


SUMMARY OF CRITICAL APPROACH OUTCOMES 


Longitudinal 

Situation 

Lateral 

Situation 

OK: 

Successful Landing 


OK: 

Successful Landing 

MA: 

Successful Missed Approach 


MA: 

Successful Missed Approach 

SL; 

Short Landing 


LO: 

Land Off Side of Runway 

HL: 

Hard Landing 


DW: 

Drag a Wing- Tip (or Engine 
Pod) During Landing 

OR: 

Overrun Runway During 
Rollout 


V : 

Land with Excessive 
Misalignment Angle 




RO: 

Run Off Side of Runway 
During Rollout 




DW ] 

V ] 

1 Drag a Wing-Tip 
and 

1 Land with Excessive Misalignment 




DW ] 

^o] 

1 Drag a Wing-Tip 
and 

1 Rtin Off Side of Rxinway 


TABLE VIII 

RELATIVE FREQUENCY OF BASIC OUTCOMES 
FOR U . S . AIR CARRIERS DURING 1 96L — 1 966 



NUMBER OF 

PERCENT OF | 


ACC IDEWTS / MC IDENTS 

ACCIDEI'ITS 

HiCIDEITrS 

Short landing 

19/18 

33 

1 2- 1 /2 

Hard landing 

12/i8 

21 

1 2- 1 /2 

Overrun runway 
( hydroplaning 
not a factor) 

1 1 /20 

19 


Drag a wing 
tip 

1/17 

2 1 

12 

Run off side 
of runway 

2/58 

1 

k 

^1-0 

Others 

12/15 

21 

9 
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From Table VIII it is evident that dragging a wing tip and running 
off the side of the runway are two relatively common incidents, but are 
only rarely serious enough to be classed as accidents. This leaves short 
landings, hard landings, and runway overruns as the more serious outcomes 
(in severity and frequency) . It is noted that these more serious outcomes 
are basically longitudinal situations. Whether this tendency will continue 
for Category II approaches is not known. 

Having digressed to define the critical approach outcomes, to separate 
them into longitudinal and lateral situations, and to look quickly at their 
relative occurrences, we can now proceed with a more detailed breakdown of 
the model that was shown in Fig. 17 . Figures I 8 and I 9 show the longi- 
tudinal and lateral approach outcome "trees,” including all of the various 
types of accidents . 

Conceptually, it is a simple step to generate "performance" trees from 
Figs. 18 and 1 9 considering the probabilities associated with each of 
the blocks shown. The resulting performance trees are shown in Figs. 20 
and 21 . It is pointed out that in Figs. 18 through 21 several implicit 
assumptions have been made. These are: 

• A missed approach is never elected when the airplane 
is within the landing window. 

• A missed approach is always elected when the airplane 
is outside the landing window. 

• Missed approaches are elected only when at the decision 
height, and they are always successful. 

Although these assumptions are not strictly valid in real life (due to 
the human judgment factor) , they can nevertheless still be made because 
they don't introduce undue bias in the computed results and because they 
involve pilot decision alone and not inherent aircraft/landing- system 
characteristics. Inclusion of such a decision element could be a future 
refinement to the basic model. 

The next step in the analysis procedure is to compute the various 
probabilities that go into the boxes in Figs. 20 and 21 . 
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Successful 
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Figiire 1 9. Lateral Approach Outcome "Tree" 





Successful 
Landing (OK) 



Probability of over-running runway, given that you 
made a successful touchdown 
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B. EQJUATIONB DEFININO APERQACH OUTCOME PROBABILITIES 


The longitudinal and lateral situations will first be considered 
separately. Then these separate results will be used to obtain the com- 
bined probabilities (as required for the cases of missed approaches and 
successful landings). 

Each of the probabilities of interest can be represented by an axea 
under a probability density distribution curve. For Gaussian distributions 
these areas are simple to determine. They depend on only three parameters: 

• The mean value of the pertinent variable 

• The rms deviation from the mean 

• The critical limits that define the accident 
(or whatever) 

The example in Fig. i:::2 shows graphically how these quantities affect the 
probability of a short landing. 

Figure then defines the mathematical expression relating the shaded 
area to a numerical probability (for Gaussian distributions). The function 
F( ) can be looked up in a table or it can be computed via an algebraic 
expression . 

Applying the relations in Fig. 25 to the exajnple shown in Fig. 22 
(assuming it to be Gaussian) gives the probability of a short landing as. 


tL = ^ 


V / 


( 1 : 


or (because the Gaussian curve is symmetric) 




PSL = (e) 

Using the above technique gives the probability of a hard landing as, 

_ Tp/ ^^TD ^'^max allowable ^ 
tHl “ \ 



I p(x jq) 


fxy: = Mean Value of x 

a-y, = RMS Deviation of x 
from Mean Value 


-*X 


TD 



Shaded area 
represents probability 
of a short landing 


Threshold of Runway (Critical limit for 

this example) 


Figure 22, Graphical Presentation of Parameters Pertinent 
to the Computation of the Prohahility of a Short Landing 


Probability of x > a is given 



Figure 25* Definition of F( ) for a Gaussian Distribution 

In a similar manner the probabilities of other touchdown accidents can 
also be expressed (assuming Gaussian distributions) . However, rather than 
pursuing this further, we will take a slightly different approach at this 
point . 

Accidents are very unlikely events and therefore accident probabilities 
are extremely small quantities. As a result, from a practical standpoint 
it is much more meaningfirl to consider the probability of a missed approach 
as a primary measure of a system^ s adequacy. This is supported by existing 
information indicating that at the Category II-B level the missed approach 
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rate can be as high as 40 percent. ThuS; for practical purposes, making 
the Category II-B window a very large proportion of the time, say 95-99 
percent, may be the most difficult part of the approach and landing sequence. 
Once this window is made, only the flare and removal of crab or wing-down 
are required to get on the ground with reasonable touchdown conditions. 
Therefore, the chance of an accident due to a short landing or an off-runway 
landing should decrease considerably in going from a Category I to a Cate- 
gory II window. However, most of the other outcome probabilities would be 
expected to be fairly similar to current operatioiia], experience. In any 
event, the considerations in the remainder of this report will be concerned 
only with determining the missed approach rate for various situations. In 
the next subsection the multidimensional approach window will be discussed. 

C . DEFIOTTION OF APPROACH WIKDOW 

The purpose of an approach window is to define the critical limits 
for continuing an approach beyond the decision height. Once these limits 
are determined it is a relatively simple matter to express the probability 
of a missed approach in a concise form. 

The approach window is defined in terms of three parameters: altitude, 

lateral displacement, and airspeed. Figure 2h shows the actual Category II 
window in space that trie airplane must fly through at the 100 ft decision 
height (in addition to maintaining airspeed within the airspeed "window”). 

The values of the variables that were selected to define the longitu- 
dinal and lateral approach windows are: 




max 


12 ft 


ye„ 


72 ft 




The d^ and y^ limits correspond to 75 microamps 1 dot) of glide slope 
error and 25 M-a (= 1 /5 dot) of localizer error at 100 ft altitude. (Full 
scale is 11 5O ,aa.) These are limits the FAA. requires for Category II 
operation in Ref. 11 . The airspeed limit was also selected on the basis 
of FAA criteria (Ref. 11 ). Note that the y^ limits correspond to the 
centerline of the airplane being v;ithin the lateral conl'ines of a typical 


1 50 ft wide runway. This means that even if the airplane is at the edge 




Figure 2^t. Category II Approach Window (at 100 ft Altitude) 
for Vertical and Lateral displacement Deviations 


of the lateral window, only a small lateral correction is required to get 
both sets of main landing gear over the runway. At this point we might 
also note that a steady 12 ft altitude error on a 2.8 deg glide path cor- 
responds to a longitudinal touchdown error of 2^6 ft, which still puts 
the airplane well within the acceptable landing zone . 

D . EQUATIONS DEFINING THE PROBABILITY 
OF A SUCCESSFUL APPROACH 

By making the assumptions of stationarity and normality (Gaussian) 
for the distributions of the random inputs (discussed in Section III), we 
can compute the means and rms deviations of the variables from their mean 
values (via the inputs and equations of motion). Then the probability 
density distribution plots, such as those shown in Fig. 25^ can be used 
to determine the overall approach success probability as follows. 

The probability of being outside the lateral limits of the approach 
window is given by the area under the tails of the lateral displacement 
distribution that falls outside the window limits- Thus, 
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Note; 

1) Longitudinal window is 
defined as Idd 

and |Au|<Aumax 

2) Lateral window is 

defined as|y«!< 

3) Shaded areas represent the 
probabilities of exceeding 
window boundaries 


Figure 25 ♦ Probability Density Distributions for Deviations from 
the Glide Slope Beam, Localizer Beam, and iJominal Velocity 


kG 


The probability of being outside of the longitudinal window is a little 
more complicated in that it involves altitude and airspeed deviations. 



Now all that remains to define the combined probability of a successful 
Category II approach is to combine the above results. Thus^ assuming that 
longitudinal and lateral deviations are independent ^ 


^MA-total ^out-long ^out-lat(”^ ^out-long) 


and 


• success 


ful approach ~ ^ "^^total 


( 19 ) 


- ( ^SA-long ) ( ^SA-lat ) 


In the next section the calculation of the various \i* s and a's will 
be discussed and some example calculations will be presented. 


47 



SECTION V 


E3CAMPLE CALCUIATIONS 


From Section IV it was found that the quantities required to compute 
the desired probahilities are the means and rms deviations from the means 
for certain of the variables . Fortunately^ the calculations of these means 
and rms deviations are quite simple to perform. The mean values of the 
variables are computed from deterministic inputs^ and the rms values 
from zero-mean Gaussian* (random) inputs^ as indicated schematically 
in Fig. 2 b. 


Deterministic Inputs 


Mean Values 

Zero - Mean 

Closed- 

Loop 

System 

RMS Values 

Random Inputs 




Figure 26. Schematic Representation of Effects 
of Deterministic and Random Inputs 

The deterministic inputs are easily handled via either a digital or 
analog computer (where the inputs are applied to the system as functions 
of time^ and time histories of the peidinent variables are recorded) . 

The value of each variable at the decision height is a number. These 
numbers are^ by definition^ the means of these quantities at the decision 
height. Clearly^ if there are no deterministic inputs^ then the mean for 
each variable is zero. 


*It is noted that Gaussian inputs applied to a linear system will 
produce Gaussian outputs for each variable. Further^ the rms value of 
each output will be directly proportional to the rms value of the input 
(for a single input'. For several independent inputs^ the rms value of 
each output is the square root of the sum of the squares of the rms values 
due to each input . 



A digital computer is used to compute the rms values for Gaussian 
inputs. The inputs are entered in the form of filtered white noise (with 
the filters determined by the input power spectra) , The computer program 
then gives the rms values as outputs. If there are no random inputs^ then 
the rms deviation from the mean for each variable is clearly zero. ThuS;, 
you get a o from zero-mean random inputs^ and a mean from deterministic 
inputs. Therefore^ for a linear system^ a combination of deterministic 
and random inputs gives a mean and a o for each variable. 

Four longitudinal examples and one lateral example will be presented 
to indicate how system comparisons can actually be made^, as well as to 
clarify the procedures involved. For this purpose a DC -8 airplane was 
selected. Of the four longitudinal systems^ three are automatic and one 
is a manually controlled flight director. All of the systems are briefly 
described in Appendix B. The three automatic longitudinal systems are 
designated as Systems A, and C. (For easy reference^ System A repre- 
sents an advanced system^ and System C represents a more conventional 
system. System B represents a system that is more sophisticated than 
but less than A.) The automatic longitudinal system examples will be 
presented first . 

For the first example, the inputs will be random gusts for Ug and Wg, 
and Gaussian glide slope beam noise. Using the closed-loop transfer func- 
tions for de (deviation from the glide slope beam) and u to "d" commands 
and gust inputs, the rms values of d^ and u for unit gust inputs (i.e., 

1 ft/sec rms) and the total beam noise are determined via Eqs . 20 through 
25- Note that new symbols are defined in these equations to simplil'y the 
writing of subsequent equations. 
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The results of tliese computations are summarized in Table IX for the 
"primary" system (A) ^ and for the other two automatic systems as well. 

It is noted that the values given in the table can be used directly to 
compare the relative merits of the three systems. Although it is not the 


TABLE IX 

SUMMARY OF LONGITUDINAL RMS VALUES FOR SEVERAL INPUTS 


SYMBOL 


LEETUITTON 


UNITS 





SYSTEM A 

SYSTEM B 



o 2 due 

to 

CVg = 

1 ft/ sec 

0.5660 

1 .052^1 

1 .1377 

ft 

Cdg 

Oj due 

to 

^Ug - 

1 ft/sec 

o.l^ino 

1 .5052 

1 .1543 

ft 

Cd 

gs 

0^ due 
noise { 

to 

10 

glide 
pa rms 

slope 

1 . 6c &5 

1 . 6155 

i • c c 

ft 

ou^ 

0^ due 

to 

Owg = 

1 ft/sec 

0-CC 

0.5527 

0.5088 

ft/ sec 

^U2 

0^ due 

to 

Oil = 

1 ft/sec 

0.501I5 

0 . 5572 

0.5498 

ft/sec 

j 

0^ due to 
noise ( 1 0 

glide 
pa ims 

slope 

) 

0,21^06 

0.211 5 

0.1714 

ft/ sec 
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intent here to compare systems on such a basis, Figs. 27 ^ 28, and 29 are 
examples of how such a comparison can be made graphically (using the num- 
bers from Table IX) . This particular comparison is made on the basis of 
random Ug, Wg, and glide slope noise applied individually to each system. 
In these three figures, parameters other than just Cu and are given 
for completeness. Several comments can be made concerning this comparison 
First, System A is seen to be considerably better than the other two sys- 
tems on the basis of (although at the expense of larger a^, and 

— but this is expected because System A has a higher bandwidth) . 
Second, on the basis of all three systems are about the same. This 
is because no considerations of speed control were made in deriving any 
of the loop closures. Thus, each system exhibits essentially the basic 
speed characteristics of the bare airframe. Further, it is noted that for 
these random type inputs applied to Systems B and C the magnitude of 
expressed as a fraction of is a lot larger than expressed as a 

fraction of %iax* This means that the probabilities of missed approaches 
for Systems B and C are almost entirely due to d-excursions , and have only 
an insignificant contribution from speed variations. This is a particu- 
larly interesting result, considering that no speed loop was closed*. Thus 
the fact that Systems B and C are considerably inferior to System A on a 
successful approach probability basis is a direct result of the larger 
d-excursions of Systems B and C, along with the fact that d-excursions 
(and not u-excursions ) produce most of the missed approaches for B and C. 

fretting back to the step by step presentation of the technique for 
computing probabilities, the ’*total" rms values of d^ and u are found from 


= 


V: 


2 2 , 2 .2 
•^di %g 


^d2^Ug ^d_ 


gs 


and^ 




2 2 2 2 2 

Q-t O,. + 0,, 0„ -f 0,-1 

U] Wg U2 Ug Ugg 


(26) 


= 7 ) 


Equations 26 and 27 can be simplified by making use of the relation 

between 0 ^ and 0 ^ given in Section III for an altitude of 100 ft. That 

g g 

is, 

= 2-590-wg 


(28) 
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Figure 27. o- Values for Random w Input of Approximately i kts 
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Figure 28 . o-Values for Random Input of Approximately 10 kts 


(deg) 





Thus, using Eq. 28 and the values in Table IX for System A, Eqs . 26 
and 27 become. 


ads(ft) = yi.62aj +2.71 

0 

(29) 

a^(ft/sec) = Jo.l^Oa^ + O.O58 

(30) 


Figure 50 is a plot of these last two relations (for System A), and shows 
that the glide slope noise contribution is negligible for greater 

o 

than about 2 ft/sec. 

At this point we could pick a value for (for example, the value 

S 

that is exceeded 1 percent of the time) and obtain values for and 
These values could then be substituted into the equations given in Sec- 
tion IV to compute the probability of a longitudinal missed approach. In 
fact, several values of could be selected and the probability of a 
successful longitudinal approach (Pg^ = 1 “ could be plotted as a func- 

tion of a-w * Such a plot is presented in Fig. 31 for System A, where it 

o 

can be seen that the probability of a successful longitudinal approach 
drops below 0-5 when exceeds about 9 ft/sec. 

Similar probabilities could be computed for alternative systems, and 
comparisons could be made. However, these probabilities are conditional 
probabilities (because they are based on the assumption of a given gust 
level) and bear no relation to the actual longitudinal missed- approach 
probability, except when the wind conditions are as assumed- Therefore, 
a further sophistication will be introduced here to account for the dis- 
tribution of gust levels . This will enable an ’’overall" probability of a 
longitudinal successful approach to be made — one that will be more mean- 
ingful in terms of long-time nationwide averages . 

This overall probability is preferred over the conditional probability 
as a performance metric for comparison of systems because it makes pos- 
sible a quantitative assessment of relative system merit. That is, a 
dollars -and- cents value of one system over another can be made if an over- 
all approach success probability is known — as opposed to only knowing 
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Figure 30. (Total) and (Total) as Function of the Wg Gust Intensity 



Figiire 5 I . ProlDability of a Successful Longitudinal Approach 
as a Function of the Wg Gust Intensity (System A) 

that one system will result in 10 percent more missed approaches than 
another when is 5 ft/sec, and 17 percent more missed approaches when 

o 

a-w is 5 ft/sec. The hey point here is that the ratio (as well as the 
S 

difference) of missed approach probabilities for two different systems is 
a function of the gust level encountered. Therefore, the relative value 
of one system over another will depend on the gust level encountered. 

This is not to say that a probability calculation made on the basis of 
a given gust level is of no use, because it is easy to conceive of a 
situation in which the decision to fly to an alternate field is made 
prior to arriving at the primary destination on the basis of a very high 
gust level at the primary destination — and its associated low probability 
of a successful approach. What we are saying is that system comparisons 
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should be made on the basis of overall probabilities in order to make 
differences between systems assessable in meanin^-ihl terms. 

The manner in which the distribution of gast levels is taken into 
consideration is seen in Eq, 51? 


MA. 


= 0.2P(m|ow =0) +0.8 


f 




)da 


Wn 


51 ) 


where 


is the overall probability of a missed 
approach 


P(MA|o^g 


is the conditional probability of a missed 
approach; and is a function of ( see 
discussion leading to Fig, j, 1 


p( is the probability density distribution of 
^ . given that clear air turbulence is 

encountered (see next paragraph / 


0.8 is the probability of encountering clear air tur- 
bulence at an altitude of 100 ft (see Section III) 


The function r>( Ovr ) is 
by differentiation. Thus^, 


determined from 


given in Section 


p(%g) 


^-(l/2)(awg/2.5)2 
5.29 ® 


(32) 


With this last equation it is now a simple matter to compute the 
overall probability of a missed approach (or a successful approach). 

It can be done by hand calculations in a few hours ^ or via a simple 
digital computer program in a few minutes. For the example case of 
System A with random '^gast’ ^gust^ glide slope noise^ the overall 
probability of a successful longitudinal approach is 0-976. However, 
by adding a Ug^g-f^ wind shear of t kt per 1 00 it from 200 ft to 1 00 ft 
the probability drops to 0*90- (Wind shear was found to be the most 
critical deterministic input.) These numbers come from Table X which 
shows a comparison of longitudinal approach success probabilities for 
two sets of inputs applied to the three example automatic systems. Also 
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table X 


LONGITULIKAL APPROACH SUCCESS PROBABILITIES FOR THE AUTOMATIC SYSTEMS 


SYSTEM 

RAIOOM Ug AND w 
GLIDE SLOPE NO 

AND 

ISE 

RANDOM Ug AND Wg, 

GLIDE SLOPE NOISE AND 
u SHEAR = !<■ KT /100 FT I'OR 
lisT 100 FT BEFORE WINDOW 

Longitudinal 
for 0^ = h ft/sec* 

OVERALL 
LONGITU- 
DINAL Pg;^ 

Longitudinal Pg_/\_ 

for ow = ^ ft/sec* 

s 

OVERALL 
LONGITU- 
DINAL Pg;^ 

A 

0.965 

0.976 

0.78 

0.90 

B 

0.58 

0.80 

0.46 

0.71 

C 

0.63 

0.83 

0.49 

0.74 


*This gust level occurs about l 8 percent of the time (Ref. 2 ; . 


included in the table is a column giving the probability of a successful 
longitudinal approach for a gust level that is equalled or exceeded about 
18 percent of the time. This column vjas presented to sho^ that the missed- 
approach rate for a moderately gusty condition is considerably greater than 
the overall average missed- approach rate. 

This brings up an interesting point. One may \vonder if there is a 
single gust level that can be considered representative, in that it gives 
the same probability as does the overall integrated gust distribution. 

For any given system the ans>/er is obviously yes . But as a practical mat- 
ter the ans-wer is unfortunately no. The single gust level that gives the 
same probability as the overall integrated gusts for System A doesn’t \iork 
for System C, and vice versa. The problem is that the relations among rms 
values, mean values, and probabilities are very nonlinear. Consequently, 
probabilities don’t scale \}ith input disturbances. As a result, every 
time the system or an input (e.g., shear) is changed, a nev; ’equivalent*' 
gust level must be found. 

A comparison of the approach success probabilities in Table X sho"ws 
that System A is considerably better than either B or C, and that C is 
slightly better than B. System A’s clear superiority vas certainly no 
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surprise. In fact^ it ^as an expected result because System A was designed 
to be an "advanced" longitudinal autopilot — optimized for the given air- 
plane dynamics. Systems B and G represent more modest autopilots (although 
the gains were optimi 2 :ed for the given airplane). However, the result 
showing System C to be better than B might be considered somewhat of a 
surprise because time traces of the "d" and "u" responses to Ug and Wg 
step inputs (given in Ref. 12) show System B to be slightly superior to C. 
On examining these time traces, it is found that the h response to these 
step gust inputs indicates a slight advantage of System C over B. The 
significance of this is explained as follows. Pirst, the d response to a 
step gust is clearly indicative of the step response of a system. Next, 
the h response to a step gust can be used to approximate the d response 
to a step gust. But the d response to a step gust is the same as the d 
response to a gust impulse. Thus the h response to a step gust is indica- 
tive of the d response to a gust impulse. The net result is that System B 
is slightly superior to C for step gust inputs, but is slightly inferior 
to C for impulsive gust inputs. Because the response to a random gust 
input more closely resembles that from an impulsive -type input than that 
from a step-t;ype input, the probability of a successful approach (with 
the random gust input) favored System C by a small airiount . However, 
because the two systems are so similar, the determination of the better 
system surely must be made on the basis of a more comprehensive comparison 
than just a random, gust input disturbance. If a step '^frust ij^put had 
been used, System B would have been found to be clearly superior to C. 
System B would also be found to be superior tor a pitch attitude bias 
input. Thus it is fully recognized that the above result favoring System 
C is only one of several criteria for judging the relative merits of com- 
peting systems. In addition to subjecting the system models to other 
inputs, such factors as reliability, maintainability, etc., must also be 
considered . 

At this point it seemis that a lesson can be learned from the above 
discussion comparing Systems B and C. This is that one must not be blinded 
by preconceived notions of what is good and what is not so good. If, in 
fact, the input that showed System C to be superior to b is encountered 
regularly, then it may be the pertinent m^etric to use in comparing the two 
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systems. In other words^ one must not judge too hastily; a system may 
have many advantages over an alternative system^ but it may be inferior 
in the one category that counts most. 

To demonstrate the applicability of the approach model to a human-pilot 
situation, the approach success probabilities for a manually controlled 
flight director approach (in the same DC-8) will be presented next. 

The mechanization of the flight director, as well as estimates for 
the pilot *s describing function and remnant, are required to define the 
manual system. These items are described in Appendix B. A key point in 
the analysis technique is that once these elements of the system are known 
(or estimated), the analytical procedures are identical to those for an 
automatic system. As a result, the manual example would involve just a 
repeat of the prior calculations with new numbers. Because this would 
not help in clarifying the use of the approach model, only the resulting 
probabilities will be presented here (to enable comparison with the auto- 
matic systems). Table XI shows the longitudinal approach success proba- 
bilities for the manually controlled flight director system with and with- 
out the estimated remnant. It is seen that the effect of the remnant is 
small (compared to the effect of the random gusts) as far as the approach 
success probability is concerned. This is a consequence of the assumption 
that the remnant scales with the displayed error (see Appendix B) . When 
the error is large there are missed approaches regardless of the remnant, 
and when the error is small the remnant contribution is also small. The 


TABLE XI 

LONGITUDINAL APPROACH SUCCESS PROBABILITIES FOR THE 
MANUALLY CONTROLLED FLIGHT DIRECTOR SYSTEM 


RANDOM Ug AND Wg AND 
GLIDE SLOPE NOISE 

RANDOM Ug AND Wg, GLIDE SLOPE 
NOISE, AND REMNANT 

LONGITUDINAL Pq^ 
FOR ft/sec* 

OVERALL 

LONGITUDINAL PsA 

LONGITUDINAL Pg^ 
FOR a-^ = 4 ft/sec* 

0 

OVERALL 

LONGITUDINAL PsA 

0.65 

0.85 

0.60 

0.81 


*This gust level occurs about l8 percent of the time (Ref. 2 ). 
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net effect is that the reimaant is not a very significant contributor to 
missed approaches. 

To complete the calculation of a combined overall missed approach proba- 
bility, a lateral example is presented next. Because the analytical tech- 
nique is again the same as that used for the longitudinal calculations, it 
is again unnecessary to repeat the computational details. Therefore, only 
the major points of the lateral example vill be presented. A brief des- 
cription of the lateral airplane dynamics and example control system is 
included in Appendix B. 

Localizer noise and random gusts for Vg and Pg were used as disturbance 
inputs. Because the rms level of the gust components at any given alti- 
tude are related as shown in Eq. 9> it is convenient for comparison pur- 
poses to continue to express the gust level by the value of cr-Wp. • For the 
lateral situation the only window constraint is the deviation from the 
center of the localizer beam. Figure "^2 shows the total rms lateral dis- 
persion (due to Vg, pg, and localizer noise) as a function of the a^g 
level. It is seen than the rms dispersions are a much smaller fraction 
of the window limit than was the case for the longitudinal situation. 

Figure 33 then shows the probability of a successful lateral approach as 
a function of the a>/rr level. For corresponding gust levels it is seen 
that the probability of a successful lateral approach is considerably 
higher than that for a successful longitudinal approach (see Fig. 31 ) • 

By integrating over all gust levels, the overall xjrohability of a success- 
ful lateral approach i. s found to be 0 - 9987 * The addition of a lateral 
wind shear of d kts per 1 00 ft from 200 ft down to 1 00 ft would give a 
lateral offset at the window of l4 ft. When the random gusts and locali- 
zer noise are superimposed on this shear the probaoility of a successful 
lateral approach is lowered to 0.9975* These probabilities are summarized 
in Table XII, which is the lateral counteipart to Table X, Finally^ when 
the lateral system is combined with each of the three longitudinal systems, 
the overall probability of a successful approach can be computed for each 
combination. These results are given in Table XIII, 

The above examples show the usefulness of the model as presented 
herein. It provides a practical means for comparing systems (on the basis 
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TABLE XII 


LATERAL APPROACH SUCCESS PROBABILITIES 


RANDOM Vg AND pg 
AND LOCALIZER NOISE 

RANDOM Vg AND pg, LOCALIZER 
NOISE AND Vg SHEAR =: 

8 KT/iOO ft for last 100 FT 
BEFORE WINDOW 

LATERAL PsA FOR 
= A FT/SEC* 

OVERALL 
LATERAL PgA 

LATERAL PsA FOR 
= 4 FT/SEC* 

& 

OVERALL 
LATERAL PsA 

0.9997 

0.9987 




*This gust level occurs about l8 percent of the time (Ref- 2). 


TABLE XIII 

COMBINED (LONGITUDINAL AND LATERAL) 
APPROACH SUCCESS PROBABILITIES 


LONGITUDINAL 

SYSTEM 

RANDOM Ug, Vg, Wg, 
AND GLIDE SLOPE AND 
LOCALIZER NOISE 

RANDOM Ug, Vg, Wg, GLIDE 
SLOPE AND LOCALIZER NOISE, 
Ug SHEAR = 4 KT/i 00 FT AND 
Vg SHEAR - 8 KT/100 FT FOR 
LAST 1 00 FT BEFORE WINDOW 

PSA for 
4 FT/SEC* 

OVERALL PsA 

PSA for 

<Jwg - ^ FT/SEC* 

OVERALL Pg^v 

A 

0.965 

0.975 

0.78 

0.90 

B 

0 .^ 

0.80 

0.46 

0.71 

C 

0.65 

0.85 

0.49 

0.74 


*This gust level occurs about 18 percent of the time (Ref- 2). 


of approach success probabilities), rather than an "academic" means (as 
•would be the case if all of the various aspects of each system were not 
combined to assess the critical performance parameter — or parameters). 

It is pertinent here to make a few comments concerning the effect of 
a nonzero mean on the probability calculations. As an example, Fig. 
shows the effect of a nonzero mean on the probability of landing off the 
side of the runway. (This example was chosen for ease of visualization.) 
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Figure 3^ . Effect of Nonzero Mean on the Probability of Landing 
Off the Side of a Runway (Represented by shaded area) 

It is seen that the area under the right-hand tail is diminished^ and the 
area under the left-hand tail is increased due to the mean being to the 
left of the center of the runway. The ijnportant point here is that the 
decrease in area on the right-hand tail is mucli less than the increase 
in the area on the left-hand tail. In fact^ it is generally true that 
when the mean moves very far from zero, the area under one tail becomes 
so large compared to the area under the other one that the smaller tail 
can be neglected. The significance of this is that vhen shear- type inputs 
(and others which give a nonzero mean) are combined with random inputs 
the probability of landing off the side of the runway, for example, becomes 
quite large compared to the probability of landing off the side of the run- 
way when only zero-mean rjaussian-type inputs are considered. (This same 
situation can exist v/ith the other types of accidents, such as long landing 
short landings, hard landings, etc., as well as for the situation at the 
approach window, v^here the parameter of interest is the probability of an 
approach success.) Thus, since zero-mean random gusts are present most of 
the time^, a relatively hxgh percentage of those occasions when shears 
(for example) are encountered will result in massed approaches (or acci- 
dents). Conversely; few missed approaches (or accidents) will occur when 
shears are not present- 
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SECTION VI 


SIM4AHY AND CONCLUSIONS 


In this report ve have confined our attention to the final approach 
phase of flight, from glide slope acquisition do™ to the Category II 
decision height of 100 ft. During this phase the primary tasks are glide 
slope and localizer tracking, i^hile the pertinent disturbance inputs are 
random and discrete gusts and beam noise (bends) . The controllers used 
in the examples included fully automatic systems, as v/ell as a manual 
flight director system. Fiirther, the example calculations included both 
the longitudinal and lateral situations. 

With the above ''components” of the approach system model it is possible 
to compute perform.ance metrics. The key feature in this report is that 
-we have derived an analytical expression for the overall approach success 
probability that is sensitive to changes in the basic airframe, the con- 
troller logic and mechanization, the displays used, types of inputs con- 
sidered, and the dimensions of the approach "window at the decision height. 
Thus, competing airframes, control systems, or displays (e.g., flight 
directors) can be compared directly. The model has a further capability 
in that the effects of changes in the approach window can easily be 
assessed. 

The numerical examples presented in Section V showed the effects of 
changes in feedback loops and equalization on the probability of a suc- 
cessful approach for a DC-8 type airplane. The example exercises showed 
quantitatively that a combination of inputs, such as wind shears and 
random gusts, leads to a significantly higher missed approach rate than 
would exist with either input by itself. The examples also showed that 
beam bends do not play a significant role in producing missed approaches. 
Further, an overwhelming majority of Category II missed approaches would 
be the result of altitude dispersions (as compared to those resulting from 
lateral position dispersions) for the advanced automatic system (System A) 
in a DC-8. 
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APPENDIX A 


ACCIDENT AND INCIDENT STATISTICS 


The air carrier terminal area accident and incident statistics for 
196^^ 1965^ and. 1966 are presented in this appendix. This information 
can be used to validate certain aspects of the approach model^ as well as 
to point out aspects of approach and landing that deserve further study. 

The necessary data for validation of the model consists of 

• Nxxmber of air carrier operations 

• Number of air carrier instrument approaches 

• Number of air carrier incidents which are takeoff- 

or landing-related 

• Number of air carrier accidents which are takeoff- 
or landing-related 

for a given interval of time. In the latter two sets of data it is also 
necessary to distinguish between landing incidents/accidents which occur 
from factors originating during instrument flight in distinction to visual 
flight- Further breakdown into type of incident/accident ^ type of aircraft^ 
etc.^ is helpful for confirming finer grain structure of the model. 

The numerical data referred to above is used for constructing the 
probabilities for occurrence of different types of incidents/accidents . 

The basis ^ or population, used for calculating these probabilities is 
one-half the number of air carrier operations, or the number of air 
carrier instrument approaches, as is appropriate. (Numbers are given 
in Tables A- I and A- II for each of three relatively recent years for 
which complete and consistent data are available.) Notice that proba- 
bilities are therefore on a per- takeoff, -landing, or -instrument approach 
landing basis, which is directly applicable for terminal operations study. 

If required, these statistics can be tied in with many of the better-known 
(but here inappropriate) statistics based upon passenger -miles, flight 
hours, or whatever, by researching appropriate conversion factors. Of 
course, it must be appreciated that terminal area incidents /accidents 
do not account for all of the incidents /accidents that occur. 
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TABLE A- I 


NUMBER OF AIR CARRIER AIRCRAFT OPERATIONS* AT FAA FACILITIES^ 


YEAR 

1964 

1965 

1966 

7,C‘T,43‘'‘ 

7,819, i44 

8,206,377 


TABLE A- II 

NUMBER OF AIR CARRIER INSTRUMENT APPROACjlES* 


YEAR 

1 964 

1965 

1966 

: 64, 195 

620,645 

66^^ , ^i55 


Tlie data summarized here concerns only the bakeoil' and landing phases 
of flight as defined by the CAB-'' The subphases of flight, as defined by 
the CAB, have been used to distinguish instrurjient fligrit incidents/accidents 
in the landing phase from the total of landing irii ‘iderits/accidents . A 
similar approach cannot be applied to the takeoff data incidents/accidents 
because no means l‘or distinguishing which arise duririg instinment flight 
exists. Furthermore^ only those air carrier incidents /accidents are con- 
sidered for which tlie aircraft involved has tv/o or more engines, weighs 
more than 15,000 lb empty, has fixed wings; and which arise because of 
piloting difficulty, deficiency, or error; ATC deficiency or error; funda- 
mental aircraft or ground facility deficiency ''suc-i as navigation facility 
failure, improper runv;ay maintenance, designed-in limitations, etc., which 
are primary for the takeoff /landing execution). Hov/ever, mechanical/ 
electrical/hy^draulic failure induced incident s/acrcidents are not counted if 


* Includes operations/instrument approaches by foreign air carriers and 
all nonfixed air carriers as well as by small and/or single-engine aircraft 
used in air carrier service . 

^ FAA Air Traffic Activity — Calendar Year 1^66 , federal Aviation Agency, 
Office of Management Services, Jan. 1 9oS / 





• The resulting take off /landing is acceptable. 

• The failure is beyond compensation by emergency 

piloting procedure^ improper maintenance 

induced gear collapse on landing, even if an 
incident or accident occurs. 

When improper implementation of emergency procedures is a heavily 
contributing cause to the incident /accident, it ^ counted. This 
particularized view sometimes leads to interpretations of causes which 
deviate from those in the official FAA/CAB report on an incident/accident . 
This is necessary in order to represent accurately those terminal opera- 
tions which are performance sensitive to the pilot- controller- display- 
aircraft- guidance system organization. 

A statistical summary of accidents for the calendar years 1964 through 
1966 constitutes Tables A-III and A- IV. 

A similar statistical summary of Incidents is also available. The 
need for incident statistics is twofold. By regarding incidents as 
"accidents wherein the deviation of the circumstances (variables) from 
nominal is not so large as to cause damage or injury, " a broader data 
base and consequently higher confidence levels in the statistics are 
gained. This also provides a means for the generation of the very lowest 
probability events which, on the basis of accident data, would have to be 
approximated by guesstimate. Secondly, incident data produce an additional 
set of conditions by which the terminal operations model may be validated 
or constructed (but not both). 

A statistical summary of terminal area aircraft incidents during 
1964 - 1966 is included here as Table A-V. The availability of these 
incident statistics broadens the small base of data on deviations from 
desired flight conditions. However, it must be pointed out that these 
new data are generally not as reliable as the data concerning accidents. 
This is due to the lack of uniformity in reporting incidents. Each air- 
line (and even different flight crews from the same airline) has its own 
point of view with regard to situations not severe enough to be considered 
accidents. In pau’ticular, less care is exercised in accurately recording 
the conditions leading to and surrounding an incident than would be the 
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TABLE A- IV 


OVERALL SM 4 ARY OF TER^IINAL ACCIDENTS 


TYPE OF OPERATION 

OCCURRENCE PROBABILITY X lo"^ 

NO. OF 
ACCIDENTS 

1954-1966 

Takeoff 

9.372 (Takeoff Acc. /Takeoff) 

1 1 

Landing (total) 

48,57 (Landing Acc , /Landing ) 

57 

(instrument flight) 

^18.67 (Inst. Land. Acc. /Inst. App. ) 

9 

Terminal 

! 

1 

28.97 (Accidents/Operation) 

68 


case for an accident. Although this does not nullify the incident 
data^ it does mean that fine grain details based on incident reports 
should not be used ”to build a story around," 

In addition to the gross numbers given in Table A-V (and the earlier 
tables) we also have a breakdown of the accidents and incidents by air- 
craft type. This data is presented in Table A-VI. The purpose for 
obtaining this information is to be able to particularize the system 
model to any given aircraft type^ and have data to help in checking 
any differences predicted between aircraft. 

It is noted that the relatively large number of landing accidents 
and incidents shown in Table A-VI for the B -727 is due^ in part^ to the 
newness of the airplane (during 196^-1966). In other words, an asymptotic 
level had not been reached on the learning curve for this airplane. That 
the DC -9 <ioes not show such a large number of accidents /incidents is 
probably due to the small number of aircraft that were in operation 
during the years of interest (compared to the number of 727s). 




TABLE A-V 
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•These data Wave been siimmarized from public records made available through the Flight Standards Division of the FAA. 

The data have been culled from the public records and classified according to the ground rules discussed in the 10 Oct- 1 966 
Quarterly Progress Report. In particular, incidents involving equipment failures only are generally excluded from consideration- 
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'Doce not meet gross weight ground rule 














APPENDIX B 


EXAMPLE AIRPLANE AND CONTROL SYSTEM CHARACTERISTICS 


This appendix contains a brief description of the longitudinal and 
lateral characteristics of the example airplane and several alternative 
approach controllers. A detailed development of these controllers (based 
on guidance^ control^ and regulation requirements) is presented in Ref. 12, 

The example airplane is a DC -8 defined by the landing approach 
configuration parameters given in Table B~I, For the purpose of demon- 
strating the use of the approach model^, three automatic longitudinal 
systems a manual longitudinal flight-director system^ and an automatic 
lateral system were used in the example calculations. 


TABLE B-I 

DC -8 PARAMETERS FOR LAM)ING APPROACH CONFIGURATION 


1 GEOMETRY AND INTERIAL 


LONGITUDINAL 


LATERAL 


PROPERTIES 


STABILITY AXES 


BODY AXES 

h 

(fM:) 

0 

Xu 

( 1 /sec) 

-0.0375 

Xv 

( 1 /sec) 

-0.0887 

M 

(-) 

0.204 

Xw 

( 1 /sec) 

0.136 . 

Xbi 

( 1 /sec) 

0 

VTo 

(ft/sec) 

228. 


(ft/sec^/rad) 

0 

Ys* 

( 1 /sec) 

0.051 


(deg) 

-2.8° 

^&T 

(ft/sec^/^) 

0.106 


( 1 /sec^) 

-1 .40 

q 

(Ih/tt^) 

6l .8 


( 1 / sec ) 

-0.283 


( 1 /sec) 

-1 .04 

s 

(«;2) 

2758. 


( 1 /sec ) 

-0.750 

Lr 

( 1 / sec ) 

0.474 

b 

(ft) 

i42.4 


(-) 

0 


( 1 /sec^) 

1.13 

c 

(ft) 

22.16 

2Se 

( ft/sec^/ rad) 

-9.25 


( 1 /sec^) 

0.159 

w 

(lb) 

180,000. 


(ft/sec^/^) 

“0.00097 

Np 

(l/sec2) 

0.368 

m 

( slugs ) 

5,580. 


( 1 /sec-ft ) 

0 


( 1 /sec) 

“O.O29 

Ix 

(slug-ft^) 

5.2 X 10^ 


( 1 /sec-ft) 

-0.00461 

Nr 

( 1 /sec) 

-0.257 

L 

(slug-ft^) 

3.8 X 10^ 


(1/ft) 

-0.00085 


( 1 / sec^) 

0 

Iz 

(slug-ft2) 

6.6 X 10^ 

Mq 

( 1 /sec) 

-0,594 

NSr 

( 1 /sec^) 

“0.368 

^z 

(slug-ft^) 

0 


( 1 /sec^) 

-0.923 




^CG 

(^ c) 

25.2 


( 1 /sec^/^) 

0.000623 




SFo 

(deg) 

50 

% 

( 1 /sec^) 

-1 .05 




0-0 

(deg) 

0.62 

Hi 

( 1 /sec) 

-0 . 1 936 
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The control equations and functions accoinplished by the three automatic 
systems are given in Table B-II, The systems are arranged from to 
in order of decreasing complexity and capability with ’"A” also standing 
for ’'advanced^' and ’’C" for ’’conventional*” All of the systems can acquire 
and maintain position on a straight line glide slope beam with well-dajuped 
path mode responses. However^ System C is not suitable for following 
higher order paths with zero steady-state error. The major distinction 
between Systems B and C is in the 9 feedback^, which is washed out at very 
low frequencies on System B and not at all on System C. The washout is 
intended to improve the Wg windproofing_, the steady- state following of 
higher order paths ^ and to remove the effects of any steady-state 0 
biases. This is achieved at the expense of a slight amount of path 
damping and bandwidth. Consequently^ the superiority of System B in 
Wp- windproofing and steady-state operations may be offset, for other 

o 

inputs, by its smaller bandwidth. 

System A is representative of the elevator axis of an advanced 
controller, typical of the forthcoming generation of low level approach 
and automatic landing systems. Appropriate feedbacks exist for all the 
functions listed for longitudinal control in Table B-II. 

All of the example systems can be improved by the addition of an 
airspeed control. However, as they now stand without such a control, 
the speed variations of the three systems are nearly identical. Further, 
this similarity would not be changed if the sajne airspeed controller were 
added to all systems. Consequently, for the sake of siinplicity, we have 
not provided airspeed control loops. 

Block diagrams corresponding to the control equations in Table B-II 
are given in Figs . B- 1 and B-2 for Systems A and C . A block diagram for 
System B would be essentially the same as that in Fig. B-2, except that 
the attitude feedback, Kq, would be replaced by the transfer function, 

K0Two^/< Two® + 1 ) • 

An analysis of a manual control situation requires several considera- 
tions not found with automatic systems. One of these is the flight director 
mechanization. A block diagram for the longitudinal flight- director system 
is given in Fig. B-3, where it is seen that features from both Systems A and B 
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CONTROL EQUATIONS AND FUNCTIONS PERFORMED BY THE E5CAMPLE AUTOMATIC LONGITUDINAL SYSTEMS 


o 


pp 


0) 

UD 

a 

f, 

!>> 

o 

c 

<y 

g. 

<u 

Vh 

'd 

o 

•H 

Ph 

0 ) 

Pi 

-p 

o 

to 


0 ) 

60 



P 





d 





o 





ft 





CO 

g 





t> 

CD 





O 










P 


P 



d 


d 



o 


o 



ft 


ft 



w 


CO 

g 



CD 


»» 

CD 





>> 



o 


o 



d 


d 



<D 


CD 





d 



ex' 


cx< 



0) 


0 



J 


I 





g 



P 




ft 

d 




o 

o 




ft 

ft 




u 

CO 




o 

d 




ft 

% 




1 (U 





•P tto 

^ q 

CD 




O d 






O 




CO 

P 


cu 



0) 


(O 


P o 

d 




ft p 

cr* 




OJ 

d> 


T 


^ a 

cO & 

ii 


p 


i <D 

1 


ft 


U 

■a 


ft 



ft 


* — ^ 


•CD 

ft 





bD 




p o 

P 




CtJ ft 

•H 




ft CQ 





P 

Q 




• *\ (D 

0 



ft 

tiO-P 




o 

P X 

ft 



p 

•rH rr} 

ft 



CO 

& 

d 

ft 



§ ^ 

•H 

P 


bD 



d 



Q ft 


ft 



ft 

ttO 




ft > >i 


d 


bO 

O ft p 

ft q ft 


(D bO 
ft d 



?H S ft 


d ft 


<<H 

0) ft ft 

(U 

O > 


o 

ft ^ 


1 O 


o 

1 ft cd 


P H 


u 

P O ft 

t 

<D i-i 


ft 

^ o cd 

ft o 

6 

ft 

O ft o 

o 

tjOft 

ft 

d 

ft 

ft 

•H 

U 

ft 

CQ 

CQ 

ft 




0 ) 

oO 

♦ 


( 1 ) 

cO 


■g 


•H 

u 

<D 

Pi O 

I ^ 

bO C 
P o 3 
O P 

!>> 

•SS 

0) 

eO cr» 

I* 

CD 


0) 

CO 


<D 

+ 




Ieh 


o 

CD 


f 


CQ 

<D 




+ 

CO 




CO 

en 


o 


JL 






•<D 

•CD 


M 




C\J 

E? 


CO 


<x> 


o 

CD 


f 


CD 

I 

•H 

-P 

*P 

"H 

•H 

<U 


'd 


ft 

CD 

-P 

CO 

bD 

d 


ft ft 
I ft 
4 ^ ft 
-P 
Q CQ 
»d 
CQ 


uoTCj.'BxnSaH pH's 
Xojcxuoo 8pnq.xq.q-V 


o' d 
o ft 

^ -p 
-P CQ 




ft d> 
ft CO 

n 

ft 

ft bO 

.S ^ 

:s 


uoxq.'Bpn^aa 
puT8 xoxquoQ qx^d 


CQ 

P 

o 


W 

ft 

o 

u 

-p 

§ 

o 


B-3 





Figure B-1 . Block Diagram of Advanced Longitudinal Approach Control System [ 
















Figure B-2. Block Diagram of Conventional longitudinal Approach Control System 
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Figure B-3. Block Diagram of Manually Controlled Longitudinal Flight Director System 









are incorporated. This particular flight director represents an advanced 
flight director from Ref. 13. 

Also required for a manual control situation are analytical descriptions 
of pilot behavior and pilot remnant. The estimated pilot describing fimc- 
tion for the flight -direct or approach task is 


= K 


/s_+^\/j^\ 

s /\s + 10j 


.-0.35s 


(B-1 


This form for the describing function is consistent with the results in 
Ref. except for the integration, which was added to represent the 
pilot's behavior in nulling the flight-director command bars (as opposed 
to being content with a steady command bar offset). The magnitude of 
time delay (0.35 sec) is slightly larger than the values given in Ref. l4 
for two basic reasons: 

• The longitudinal flight-director tracking is not 
a full-attention task as was the case with the 
Ref, experiments, 

• Deviation of the manipulator from the ideal 
(heavily spring-restrained) will degrade the 
pilot's capabilities somewhat. 

An estimate of the pilot remnant was based on data taken in DC-8 
simulator experiments described in Ref. TO. It was assumed that the 
magnitude of the remnant scales with the displayed error (i.e.^ flight- 
director commands) such that the ratio of the correlated error power to 
the total error power is 0.8. This value of 0.8 was obtained from the 
reduced simulator data. No residual remnant was considered. The 
spectral shape of the estimated remnant is given by 


®nne 


a 


2 

e 


K 

s+3 


2 


(B-2) 
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Knowing the mechanization of the flight director^ and having estimates 
for the human pilot describing function and remnant allows the manual 
system to be analyzed in the same manner as automatic systems. 

Table B-III summarizes the numerical values of the gains^ time 
constants^ etc.^ for each of the longitudinal example control systems. 

TABLE B-III 

SUinvIARY OF NUMERICAL VALUES DEFILING THE 
LONGITUDINAL EXAMPLE COmOL SYGTEl^lS 


SYSTEM 

A 

B 

c 

FLIGHT 

DIRECTOR 

UNITS 


-2.0 

-3.652 

-3.652 

-0.615 



-2.0 

0 

0 

-0.615 



^.0256 

0 

0 

-0.00678 

rad 

ft/sec 

% 

-<'-.00R6T 

-0.005 i 4 

-0. 0051^1 

-0.001356 

rad 

TT 

% 

-0.000768 

0 

0 

0 

rad 

ft-sec 

1 

T 

-^wo 

0.7 

0.08 

0 

0.7 

-1 

sec 

1 

Ta 

1 5.0 

15.0 

15.0 

18.0 

-1 

sec 

1 

Tf 

2,0 

2.0 

2.0 

2,0 

-1 

sec 




— 

— 

— 

— 


The automatic lateral control system used in the example calculations 
is defined in the block diagram of Fig. B-t. As with the longitudinal 
systems^, guidance^, control^ and regulation requii-ements were the primary 
considerations that went into the development of the lateral system. 
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i^;^ure , Block I'iagrajn of :..ateral .Approach Control System 










